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Abstract 

A numerical method is developed to calculate the transition tempera- 
ture of double or multi-layers consisting of films of super- and normal con- 
ductors. The approach is based on a dynamic interpretation of Gorkov's 
linear gap equation and is very flexible. The mean free path I of the differ- 
ent metals, transmission through the interface, ratio of specular reflection 
to diffusive scattering at the surfaces, and fraction of diffusive scattering 
at the interface can be included. Furthermore it is possible to vary the 
mean free path and the BCS interaction NV in the vicinity of the inter- 
face. The numerical results show that the normalized initial slope of an SN 
double layer is independent of almost all film parameters except the ratio 
of the density of states, (d s /T s ) \dT s /dd n \ = F 3n (N n /N a ). There are only 
very few experimental investigations of this initial slope and they consist 
of Pb/Nn double layers (Nn stands for a normal metal). Surprisingly the 
coefficient r sn in these experiments is of the order or less than 2 while the 
(weak coupling) theory predicts a value of about 4.5. This discrepancy 
has not been recognized in the past. The autor suggests that it is due to 
strong coupling behavior of Pb in the double layers. The strong coupling 
gap equation is evaluated in the thin film limit and yields the value of 1.6 
for r 3n . This agrees much better with the few experimental results that 
are available. 

PACS: 74.45.+r, 74.62.-c, 74.20.Fg 
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1 Introduction 



The transition temperature of a thin superconducting film in contact with a 
normal metal is reduced. This is known as the superconducting proximity ef- 
fect (SPE). The double layer SN or a multi layer (SN)„ can consist (i) of a 
superconductor S and normal conductor N or (ii) of two superconductors with 
different transition temperatures (the one with the lower T c is generally denoted 
as N). Its systematic experimental investigation started in 1960's by the Hilsch 
group in Goettingen PP, El and stimulated a number of further experimental 
investigations |3], @], IE El- For the dirty case (mean free path of the con- 
duction electrons is much smaller than the coherence length) Werthamer 
derived a set of implicit equations for the transition temperature of double lay- 
ers consisting of two superconductors. After some modification according to de 
Gennes' boundary condition 2} between the superconductors, the Werthcimcr 
theory described the experimental results for double layers of two superconduc- 
tors quite well (see for example |S], |8], E|). The Wertheimer theory is restricted 
to short mean free path, (using the diffusion approximation) and uses what is 
now called the single mode approximation (the gap function is approaximated 
by a cos (fc s z)-dependence). Theoretical results for the clean case where the 
mean free path / is larger than the BCS coherence length £bcs are more diffi- 
cult and the case where 1,^bgs and the film thicknesses are of the same order 
of magnitude are much more challenging. 

In recent years the superconducting proximity effect has experienced a re- 
newed interest. A large number of papers studied the SPE theoretically [TU] . 
EE EE EH, EE Hi, EE EE EE EE EE EH, EE ESI and exper- 
imentally El, ESI, ESI, EE EE EE EE EE ESI, EE E2 particularly 

during the last 10 years. The studies have been extended to double layers of a 
superconductor and a ferromagnet (SF) [33], [3J], [55], 

Recently our group revisited the superconducting proximity effect using it as 
an experimental tool [301 • One interesting information the SPE provides is 
the transparency of the interface between the two metal films for the conduction 
electrons. The reduction of T c in the superconducting component of the SN 
double layer depends on the rate at which electrons can cross the interface 
between S and N. This interface transparency is of interest in a number of other 
disciplines and applications in solid state physics. 

When our group tried to compare the experimental results for the transi- 
tion temperature with theoretical predictions we found that only a few recent 
theoretical investigations calculated the transition temperature of SN double 
layers EE EE ESI- These papers considered the extreme cases, either the 
clean limit for infinitely large mean free path or the dirty limit 36: where 
the mean free path is much shorter than the BCS coherence length. Reference 
pTfT) considered superconductor- ferromagnet double layers in the " dirty limit" . 
It includes the case of an SN double layer by setting the exchange energy in the 
ferromagnet equal to zero. A multi-mode expansion of the order parameter is 
used in the superconductor. This yields a complex set of equations which con- 
tain the transition temperature implicitely. Their single mode approximation is 
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similar to Wertheimers result. 

Since our experiments used films with short and large mean free paths the 
author prefered to develop a numerical procedure which is capable of calculating 
the transition temperature of arbitrary sequences of superconductors and normal 
conductors in a wide range of the mean free path. This calculation uses a 
simple interpretation of the gap equation which was stimulated by deGennes 
work Below I will sketch the (simple) numerical procedure. In chapter II 
the theoretical background is reviewed and the numerical procedure discussed 
in detail. In chapter III some of the numerical results are presented. In the 
discussion of chapter IV I will point out a discrepancy between all experiments 
I am aware of which study the change of T c of a superconducting film when 
covered with a thin normal conducting film, i.e. the normalized intial slope 
d dT 

— — -, where d s and T s are the thickness and transition tmperature of the 
T s dd n 

dT 

superconductor and — - is the initial slope of the T c -reduction for zero thickness 
dd n 

d n of the normal conductor. 
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2 Theoretical Background 



2.1 The linear gap equation 

The superconducting phase transition in zero magnetic held is generally of sec- 
ond order. Therefore, close to transition temperature T c of the double layer, 
the gap function A (r) , which is the order parameter of the phase transition, 
is small and only terms linear in the gap function contribute. This linear gap 
equation, hrst formulated by Gorkov |37| was rewritten by deGennes jS] as 



Here A (r) is the gap function at the position r, w n = (2n + 1) irksT/h are 
the Matsubara frequencies, V (r) is the effective electron-electron interaction 
at the position r. The sum is limited to the range \ui n \ < where fie is 
the Debye temperature. This corresponds to a sum over n from — n c to +n c 
where n c = Qd/ (27rT) = £Idt~t, where Qjj and tto are the Debye temperature 
and frequency and tt — fij (27t/cbT). The function H LOn (r, r') is the product 
of two Green functions G UJn (r, r') and G* (r, r') which represent a Cooperon. 
Since the Green function G Wn (r, r') represents the amplitude of an electron 
traveling (at finite temperarture) from r' to r the product G Wn (r, r') G* (r, r') 
describes the amplitude of a Cooperon traveling from r' to r. Since the two 
single-particle Green functions are conjugate complex to each other, the product 
of their amplitudes is proportional to the probability of a single electron to 
travel from r' to r. If one interprets in equation (JIJ G lVn (r, r') G* (r, r') as 
the propagation of single electrons then one has an equivalent problem and its 
solution is also the solution of the gap equation. In the following the solution 
of the equivalent problem will be persued. 

From the properties of the Green functions G Un (r, r') (see appendix 6.1) it 
follows that H Un (r, r') is the electron density if one injects continuously elec- 
trons with a rate N/tt at the point r', while their density decays along the path 
as exp (—2 \ui n \ s/vp) where s is the distance traveled (not the distance from r') 
and N is the BCS density of states. 

The right side of equation (JTJ d 3 r'H LJn (r, r') A (r') (excluding £"=-n ) 
yields the density of electrons at the position r when one injects constantly 
NA(r')d 3 r'dt'/TT electrons in the incremental volume d 3 r' at the position r' 
per time interval dt' , which decay during their propagation with the decay rate 
of 2 \u> n \ (tt = h/ finksT)). (NA (r') d 3 r' represents a (dimensionless) number 
of electrons and the rate of injected electrons per volume is NA (r') /tt). These 
electrons propagate with their Fermi velocity from r' to r, either directly or dif- 
fusively. Their density decays along the path as exp (—2 \uj n \ t' A ) where t' A is the 
time since the departure from r'. At the position r the surviving density of all 
arriving electrons is integrated over J d 3 r' J_ dt' . When summed over ui n and 




|w n |<n D 



(1) 



H Un (r, r') = k B TG^ (r, r') G* Un (r, r') 
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multiplied with the attractive electron interaction V (r) one has to recover the 
original A (r) . 

For further treatment we define the propagation density p(vf', r,0; r', t' < 0). 
If an electron with Fermi velocity vp is introduced at the time t' < at the 
position r' then p (vp; r,0; r', t') describes the probability to find the electron at 
the time at the position r. With this definition we can express H UJn (r, r') 

r° dt' 
H Un (r, r') = N (r') / p (v F ; r,0; r', t 1 ) exp (-2 |w n | \t'\) — 

where l/r T = 2-Kk B T/H. 

The sum over u> n in Q applies only to the exponential decay functions 
exp (—2 \uj n \ \t'\) and yields the time function i]t (f) 

,.,s / I I ,.m 1 - exp (-2 (n D T T + 1) \t'\ /t t ) , , 
7r(*)= L exp I) = Sinh( | t ,| /TT) ( 2 ) 

|u„|<!1d 



(r2r;=Debye frequency). Then one can express the gap equation as 

A (r) = V (r) / d 3 r'N (r') / — y T (*') P («*■; r,0; r', i') A (r') 

J J-oo T T 



(3) 



It is obvious that the superconducting properties of the system occur only in 
the effective interaction V (r') and the decay function rjr (f). Of course, A (r) 
is the superconducting pair amplitude but in equation J3J it is just the eigen 
vector of the integral kernel. The self-consistancy condition requires that this 
eigen value is equal to one. 

This interpretation of the gap equation yields a natural extension to a time 
dependent pair amplitude or gap function. One obtains 



A (r, t) = V (r) / dViV (r') / — m (t') p (v F ; r,0; r', t') A (r', f 

J J-oo T T 



(4) 



From this equation one can derive a time dependent Ginsburg-Landau equation 

m 

For a homogeneous superconductor one has a constant energy gap. In this 
case one can perform the integral over d 3 r', using J d 3 rp (vf; r,0; r', t') — 1 and 
dividing by A 

1 = (VN) S f — VT (t') (5) 

which yields 

1 n " 1 
J_ = V 1 

n + \ 

n— ^ 

The condition JSJ) is used to determine the BCS coupling strength {NV) S .It 
has the advantage that it is not restricted to integer values of n c = Qdtt- 
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2.2 The gap equation for double and multi-layers 

Now we can apply the gap equation © to the proximity effect. The direc- 
tion z is chosen perpendicular to the multi-layer and the films are treated as 
homogeneous in the x-y plane. If there is no magnetic field then the gap de- 
pends only on the z direction. Therefore one can perform the integration over 
/ dx'dy , p(v F ;r,0;r',t')=-p(z,0;z',t'). 

Now the function p (z,0;z', t') describes the density at the time t = and the 
position z integrated over the x and y directions. For the numerical procedure 
it is more convenient to shift the time integration from the range (— oo,0) to 
the range (0, oo). 

dz'N(z') — r lT (t)p(z,t;z / ,0)A(z') 

Jo T T 

The multi-layer will be divided into small sheets parallel to the film surfaces. 
The layers are indexed by v and possess a thickness A„. 



N 



1 



Fig.l: A double layer of a superconductor S and a normal 
conductor N. The two films are split in thin paralllel layers v 
with the z position z v 



In the present paper we determine the gap-function A (z) at the transition 
temperature of an SN (superconductor/normal metal) double layer. We proceed 
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with the following steps which are demonstrated by Fig. 2. 
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Fig. 2: The double layer is sliced into sheets of 
thickness X s ,X n parallel to the film planes. 

• The superconductor is divided into Z 8 layers of thickness \ s where \ s = 
d s /Z s (d s is the thickness of the superconducting film). 

• The BCS interaction V s for the superconductor (s) is fitted, using the den- 
sity of states N s and the Debye temperature 6z)(appendix 6.2.1). 

• The time interval = 2X s /vp^ s is the time step of the numerical calcula- 
tion (vf, s is the Fermi velocity of the superconductor) (appendix 6.3). 

• For the normal conductor (superconductor with lower T c ) the same time 
step is used by dividing its thickness in layers of thickness X n — VF,nTd/2. 

• An initial gap function A„ = A (z„) is introduced. Each cell is occupied 
at the time t' — with 0„ (0) = N (z u ) A„A (z„) electrons. (N {z w ) is the 
local density of states, i.e. equal to N s in the superconductor) (appendix 
6.2.2). 

• A procedure for diffusive and ballistic propagation of electrons in the dif- 
ferent films is derived (appendix 6.3). 

• The maximal transmission of an electron through the interface in each 
direction is calculated. It can be scaled down to include a barrier at the 
interface (appendix 6.4). 

• The density O v (m) is calculated in descrete steps for the time t' = rriTd- 
(appendix 6.3). 
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• Due to thermal dephasing this density is, at each step, multiplied with the 
time factor t\t (jnTd) ■ 

• The sum J2 m ®v {m) r\T is formed, multiplied with (t,i/tt)/X v and, 
in the superconductor(s), multiplied with V s , the attractive electron-electron 
interaction. 

• The resulting function A„ is the input A^ for the next iteration. 

• Since the eigen value has to be 1 the ratio r = A (z v ) / ^ v A (z v ) is 
calculated. If r > 1 (r < 1) one increases (lowers) the temperature. 

• The interation process is completed when initial and final Aj, agree with a 
relative accuracy of 1CP 5 . This is generally achieved after a few iterations. 

All the step of the numrical procedure are described in details in the ap- 
pendix. 
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3 Results 



There are numerous parameters in the superconducting proximity effect: the 
coherence lengths £ Si „ = vpTx (for the superconductor this is the BCS £bcs if 
one uses the transition temperature in tt = h/ (27rfcsT), the mean free path l SjTl 
and the film thickness d S:1l for each film. In addition one has the interface and 
the boundaries. Any barrier between the two metals will reduce the transfer 
through the interface. Furthermore one can have additional scattering at the 
interface between the two films due to a mismatch of the two lattices. The two 
surfaces with the vacuum can reflect or scatter the incident electrons or anything 
in between. All these scattering parameters influcence the propagation of the 
electrons and therefore the transition temperature of the double layer. In the 
numerical calculation all these parameters can be included if they are known or 
used as fit parameters. 

3.1 Transition temperature 

In the majority of experiments the onset of superconductivity is measured for 
a double layer of a thick normal conducting film which is covered with a super- 
conducting film of increasing thickness. Therefore the first plotted numerical 
result represents a double layer of an infinitely thick normal conductor which is 
covered with a superconductor of increasing thickness. Among the large number 
of possible parameters the following choice is made: (i) the electronic proper- 
ties (iV Sj „,VFs,n) of the normal metal and the superconductor are identical, (ii) 
the mean free path of the normal conductor is infinite, (iii) the thickness of 
the normal conductor is infinite, (iv) the interface is perfectly transparent, (v) 
for the mean free path of the superconductor the following values are chosen: 
l s = oo, £o> Co/10, Co/100. The results are shown in Fig. 3a. The parameter a is 
defined asa = l s /S,o- The curves of the transition temperature versus thickness 
of the superconductor show the typical behavior; they approach T s for large d s 
and show a steep decline at a critical thickness d cr . The value of the critical 
thickness decreases strongly with decreasing mean free path l s of the supercon- 
ductor. It might be surprising that even a mean free path l s = Co shifts the 
T c -d s curve considerably to smaller thicknesses. For the smallest mean free path 
of l s = C/100 the critical thickness is about d cr « 0.19£o- In Fig. 3b T c is plot- 
ted versus the reduced thickness d s /d cr - The points lie almost on an universal 
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curve, particularly those for smaller l s . 
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Fig. 3: The reduced transition temperature T c /T s for 

an NS double layer where d n — oo, l n = oo as a 
function of thickness d s /£ (£=BCS coherence of S) a) 
For different mean free paths l n of S, the parameter 
a = | a /£. 
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Fig. 3b: The same plot as a) with the S thickness 
scaled with the critical thickness d c r. 



3.2 Pair amplitude 

In the next step the actual dependence of the gap function on position is of 
interest. In Fig.4a-d this gap function A (z r ) is plotted as a function of 2' = z/d s . 
We choose double layers where T c lies in the steep decline of the T c curves in 
Fig. 3 at about T c /T s « 0.3. Fig. 4a shows A (V) / Ao for the superconductor with 
ls/£,o — 0.01. (Since the amplitude of A (z r ) approaches zero at the transition 
temperature the value Ao at the maximum is of no phyisical significance) . Since 
the gap function has a horizontal slope at the free surface a comparison with a 
cosine function cos (p (1 — z')) is useful. The resulting fits are shown in Fig. 4a- 
d. Fig. 4a for Z s /£o = 0.01 shows an almost perfect quarter of a cosine function 
with p = 1.57 which is as close to w/2 as it can be. For Z s /£o = 0.1 the shape 
of the gap curve is still quite close to a cosine function but the factor has the 
value p = 1.46. For l s /£o = 1 the shape of the gap curve shows already clear 
deviations from a cosine curve and the coefficient is p » 1.25. Finally in the 
clean limit the gap function curves stronger for small z' than the cosine curve 
and the coefficient is p « 1.05 for the shown fit. This behavior is interesting 
because in a number of theoretical papers the gap function is expanded into a 
series (see for example ref. [SHI where a series consisting of cos (fio (d s — 
and cosh ((Cl m (z — d s ) /£<*)) is used, £d = v 7 ^? is the superconducting diffusion 
length and £lo,£l m are coefficients defined in that work). 



10 





Fig.4a-d: The gap function A (z) is plotted versus the position z/d s in the 
superconductor for NS double layers. Each drawing corresponds to one 
of the curves in Fig. 3a. d s is close to the critical thickness d cr . The ratio 
l s /£o is noted in the figures. 
The simple form of the gap function in the case of l s /£,o — 0.01 makes it very 
obvious why the very disordered superconductors (often discriminatingly called 
dirty superconductors) are much easier to describe. This becomes still more 
obvious if one compares the shape of the gap function at different T c /T s values 
(which means, of course, using different thicknesses of the superconductor). In 
Fig. 5a the (normalized) gap functions for T c /T s values of about 0.3 and 0.9 are 
shown as a function of z/d s . They lie perfectly on the same quarter of a cosine 
function. This is very different for the clean limit where the shape depends 
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strongly on the temperature. 
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Fig. 5a, b: The gap function A (z) is plotted 
versus the position z/d s for two NS double 
layers, each at two different d s (resulting in 
different T c of about 0.3 and 0.9 T s .) a) dirty 
limit Z s /£=0.01. 




b) clean limit l s = oo. 



3.3 Dirty limit 

Since in the dirty limit the gap function approached such a simple form for a 
superconductor in contact with an infinite clean normal metal it is worth to 
check the situation when both metals are dirty. This is the case which most 
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theoretical papers investigate. 
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Fig.6a,b: The function A/ (NV) for an SiS 2 double 
layer, T c l = T s 2/2 (details in text and figures). Both 
films are in the dirty limit, text and figures). Both 
films arc in the dirty limit. In a) the transmission 
coefficient from S2 to Si is t = 1.0 whil in b) t = 0.8. 
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In Fig. 6a A/ (NV) is plotted for a double layer of two superconductors 
with different transition temperatures of T n = 5K and T s = 10K. In ad- 
dition the density of states for superconductor N (with the lower transition 
temperature) is larger by a factor of 1.5 than for superconductor S. Therefore 
the superconducting coherence lengths £os = frvp s / (2irkBT s ) = lOOnm and 
£on = hvF n / (2'KksTn) = 300nm are different (The additional factor of two 
stems from the ratio of the transition temperatures). The difference in the den- 
sity of states and the Fermi velocity of the two metals yields a ratio of the two 
transmission coefficients at the interface Tjv->s /Ts->n — -AAA. 

The thickness of each film is d n = d s — 20nm. In Fig. 6a the mean free 
paths are chosen in both films to be l s = l n = lnm. For the corresponding 
superconducting diffusion lengths £ds>£<in one finds £ds = V^Osh = 10nm and 
£,dn — V^OrJn = 17. 3jwti. According to de Gennes the function A/ (NV) should 
be continuous at the interface. As can be easily recognized from the plot in 
Fig. 6a this condition is well fulfilled. Werthamer |7] expressed the z dependence 
of the gap function A (z) in the two superconductors as 

cosh (k n (d n + z)) in the superconductor with T n in the range —d n < z < 
cos (k s (d s — z)) in the superconductor with T s in the range < z < d s 

Fig. 6a shows for z < a fit to the function a cosh (k n (d n -I- z)) and for z > to 
the function cos (k s (d s — z)). The fitted curves lie within the trace of the points. 
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The fitted values for the parameters are k s = 0.0506nm _1 , k n = 0.0468nm _1 
and a = 0.351. This yields for the value of A/ (NV) on the left and the right 
side of the interface: 0.516 and 0.530. The corresponding slopes on the left and 
right side of the interface are: 4.29 x 10~ 2 and 1.77 x 10~ 2 . According to de 
Gennes the derivative (D/V) dA/dz should be continuous at the interface for 
the dirty limit. Using the input data of the two superconductors D Syn and V s ^ n 
one obtains for the ratio of the slopes 2. 61. The simulated A (z) yields a slope 
ratio at the interface of 2.42. So the de Gennes condition is verified with an 
accuracy of about 10%. 

In a second simulation the transmission through the interface is reduced by 
a factor 2. It is quite remarkable that this changes the transition temperature 
only from T c = 7.6K by just O.IK to 7.7K. In Fig.6b the function A/ (NV) 
is plotted for the double layer as a function of z. One recognizes that now 
A/ (NV) is no longer continuous at the interface. The functional form in N 
and S can still be well fitted by a hyperbolic cosine and a cosine function. (The 
fitted curves lie within the width of the numercal points). 

3.4 Initial slope 

When one condenses the normal metal on top of the superconductor then the 
transition temperature of the double layer decreases. Here the focus is on the 
question how the intial slope at d n — depends on various parameters, such as 
the mean free path in the superconductor and the normal conductor and the 
transparency of the interface. 

The dependence of the initial slope on the mean free path is shown in 
Fig. 7a, b. In both figures the thickness of the superconductor is equal to the 
BCS coherence length £o-The transition temperature T c /T s is plotted versus 
the thickness of the normal conductor. In Fig. 7a the mean free paths in both 
films are equal and vary between l s — l n — £o/10, £o and 10 3 £o- In Fig. 7b 
four different combinations of (l s ,l n ) are chosen. From the top to the bot- 
tom (l s ,l n ) is equal to (£ /100,£o/100), (&/100, 10%) , (l0 3 £ ,£o/100) and 
(l0 3 £o, 10 3 £o)- For all curves the initial slope is identical. (In all the numerical 
calculations which were discussed so far the two density of states are assumed 
equal N s = N n .) 
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1.0- 




of <i„/£. The thickness of the superconductor is 

equal to the BCS coherence length £. (N s = N n ). b) Different combinations of the mean free 
a) The mean free paths l s = l n are parameters. paths are used as parameters. 



In Fig. 8 the dependence of the initial slope on the thickness of the super- 
conducting first layer is tested. The graph shows the dependence of T c /T s for a 
small range of the thickness d n of the normal conductor to emphasize the initial 
range. 
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Fig. 8: T c for SN double layers as a function of d n /£. 
The parameter d s is the thickness of the 
superconductor. (N s — N n , l s =l n = £/10). 

In table I the normalized initial slope is collected. (The numerical points had 
to fitted with a polynomial to extract the slope from the numerical results). 
Up to a thickness of d s = £ the S sn is constant within about ±1%. For larger 
d s it decreases slightly. But since the value of dT c /dd n becomes quite small 
this thickness range is not well suited for the experimental determination of 
the slope. The main result is that the normalized initial slope is essentially 
independent of the thickness of the superconductor. 
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Table I: Normalized initial slope for different 
thicknesses d s of the superconductor 

Finally Fig. 9 shows that the initial slope does not depend on the transmission 
through the interface. In this calculation the density of states in both metals 
is chosen to be equal N s — N n and the mean free paths are l s = l n = £/10. 
The transmission coefficient is varied between 0.2 and 1.0. The resulting T c -d n 
curves show the same initial slope. 
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Fig. 9: T c for an SN double layers as a function of 
d n /£. The parameter t is the transparency of the 



of 

i 

interface. {N s = N n , d s = £, l s = l n = C/10). 
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4 Discussion 



The intention of this paper was to develop a convenient numerical procedure for 
the superconducting proximity effect so that graduate students could instantly 
compare their experimental results with the theory. One importat result of this 
investigation is the fact that the (normalized) initial slope of an SN double layer 
is independent of most film parameters except the density of states ratio and 
the effective BCS interaction. 



ds_ 

Tco 



dT c , 



dd n 



= r,„f (6) 



In the case of a weak coupling superconductor T sn is given by the Cooper 
limit, i.e. T sn — 1/ (NV) , the inverse of the BCS interaction. If the Debye 
temperature is not several orders of magnitude larger than T s then one has to 
determine T sn in equation © numerically. Using T s — 7.2K for Pb then the 
prefactor is about 4.5. (This is actually the value for a wide range of the Debye 
temperature between 100 A and 300A"). 

Recently our group investigated the proximity effect between the supercon- 
ductor Pb and several alkali metals ■ It was a great surprise that the exper- 
imental initial slope of these SN double layers could not be explained with the 
density of states from the literature. Instead the experimental (d s /T s ) (dT c /dT) 
was too small by more than a factor of two. Table II gives some of the data 
of the SN double layers. (The thickness of the normal metal was the smallest 
thickness in a full curve.) 

metals d s (nm) d„ (nm) S s „| cxp N„/N s ratio 

Pb/K 12.9 2.04 0.423 0.223 1.90 

Pb/Na 13.9 2.18 0.546 0.300 1.82 

Pb/Ag 17.9 2.10 0.625 0.335 1.86 



Table II: The normalized initial slope of SN double layers with Pb as 
superconductor and different normal metals. The columns 2-6 give the 
thickness of the superconductor, the normal conductor, the experimental 
initial slope, the ratio of the density of states and the ratio S s „| cxp / (N n /N s ). 

We searched the literature for other measurements of SN double layers and 
their initial slope. It turned out that there are very few measurements of SN 
layers. (At this stage we excluded transition metals because they show two- 
band superconductivity and it is not obvious how the different superconducting 
bands couple to the normal conductor). There were essentially two groups of 
publications which had measured SN double layers which contained informa- 
tion about the intial slope. The first group of papers was by Hilsch et al.pQ, [2J 
who investigated quench condensed PbCu layers. The second work was by Min- 
nigerode [Hj who also investigated PbCu layers but prepared the layers at room 
temperature. Particularly the second paper gives detailed tables of thicknesses 
of the two components and transition temperatures. The results of these papers 



18 



are collected in table III. The first column gives the components of the SN dou- 
ble layer, the second and third columns the thicknesses of the superconductor 
and normal metal. The fourth column contains the experimental normalized 
slope and the fifth gives the ratio N n /N 8 . The last column contains the ratio 
of the experimental (normalized) slope to the density of states ratio. Again the 
experimental normalized slopes are much smaller than the theory predicts. 



metals 


d s (nm) 


d„ (nm) 


Ssn |exp 


Nn/N g 


ratio 


Pb/Cu 1 


10.0 


10.0 


0.542 


0.448 


1.21 


Pb/Cu 1 


15.0 


10.0 


0.500 


0.448 


1.12 


Pb/Cu 2 


22.9 


3.30 


1.110 


0.448 


2.48 


Pb/Cu 2 


24.3 


3.90 


0.935 


0.448 


2.08 


Pb/Cu 2 
Pb/Cu 2 


32.9 


4.10 


0.883 


0.448 


1.97 


27.2 


13.10 


0.683 


0.448 


(1.52) 


Pb/Cu 2 


28,0 


26.40 


0.842 


0.448 


(1.88) 


Pb/Cu 2 


33.4 


17.70 


0.608 


0.448 


(1.36) 



Table III: The normalized initial slope of PbCu double layers, data 1 from 
ref. are quench condensed and data 2 from ref. [BJ are 
condensed at room temperature. The columns 2-6 give the thickness of 
the superconductor, the normal conductor, the experimental initial slope, 
the ratio of the density of states and the ratio S sra | e xp/ (N n /N s ). 

It is rather amazing that this fundamental discrepancy between experiment 
and theory has not been realized. What is the reason for this disagreement? 
The authors best guess at the present time is that the use of the weak coupling 
theory of superconductitivity is not adaquate for the double layers containing the 
supercondutor Pb. The superconductor Pb is a convenient component of an SN 
double layer because it has a rather large T c and is easy to condense. However, 
Pb is a strong coupling superconductor. The Fermi sphere of free electrons is 
modified by the electron-phonon interaction. An electron k which lies below the 
Fermi energy within the Debye energy ks&D c & n emit a virtual phonon and 
perform a transition into a state k' above the Fermi energy. As a consequence 
the states below kp have an occupation less than 1 and the states above kp have 
an occupation larger than 0, even at zero temperature. At kp the occupation 
does not jump from one to zero, but has a smaller step of 1/ (1 + A), where 
(1 + A) is the electron-phonon enhancement factor. The parameter A can be 
calculated with the Eliashbcrg function a 2 F (ft) which desribes the probability 
of an electron to emit or absorb a phonon of energy Ml. (A nice treatment 
of the electron-phonon physics can by found in the book by Grimvall 39 ). If 
one moves now an electron from the Fermi energy to a k-value above the Fermi 
energy then it changes the occupation of the k-state only by (1 + A) 1 . The 
change in kinetic energy is therefore reduced by the same factor; the energy 
states lie closer together by the factor (1 + A). As a consequence the density of 
states is enhanced by (1 + A) and the Fermi velocity reduced by 1/ (1 + A). This 
enhancement is restricted to a small region in the vicinity of the Fermi energy. 
The new density of states N* = N (1 + A) = N = m (1 + A) k F / (2ir 2 h 2 ) and 
Fermi velocity v* F — vp / (1 + A) = Kkp jra (1 + A) are called the renormalized 
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density of states and Fermi velocity and marked with a star. The Fermi wave 
number kp does not change and both changes in N and vp can be expressed by 
a change of the mass m* = m (1 + A) (therefore the name mass enhancement) . 

Many properties of the electron gas and the strong coupling superconductor 
can be reasonably well calculated in a renormalized strong coupling description. 
For the superconductivity this means that one applies the weak coupling BCS 
theory using renormalized density of states and Fermi velocity. The error for 
Pb is generally less than 20% as the ratio of 2Ao/fcsT c = 4.2 instead of 3.5 
demonstrates. (For other parameters see @U]). 

In strong coupling superconductors the effective interaction (NV) S is re- 
placed by (A — n*) where A describes the strength of the electron-phonon in- 
teraction and fi* is the Coulomb pseudo-potential. In the literature one find 
A-values for Pb which vary between 0.8 and 1.6 39 . Depending how much the 
strong coupling theory modifies the expression for the initial slope it could alter 
dramatically the theoretical prediction for the initial slope. In the appendix we 
show that a renormalized strong coupling theory would yield an initial slope 
of T sn = (A — //*)/ (1 + A). This appears to alter the theoretical value for the 
initial slope. The difficulty is that the weak coupling treatment just a replaces 
(NV) S by (A — /j,*) I (1 + A). This means that for the transition temperature 
one obtains the equivalent condition 



This does not alter the situation and yields, as observed by Morel and Anderson 
|41| . a too small value of 0.4 for A. 

An obvious proposal would be to solve the superconducting proximity effect 
for strong coupling superconductors. This means to develop and solve a series 
equations for the energy and position dependent gap function A(r,w;) which 
has the equations and J7J) as limiting cases. This is a very demanding job 
which goes beyond the scope of the present paper and has to be left for future 
investigations. However, we can check whether this extension is a promissing 
one. For this purpose we consider the analogy to the Cooper limit for strong 
superconductors . 

Let us consider a double layer composed of (very) thin films of a strong 
coupling superconductor and a normal metal. The super- and normal conduc- 
tor have A values A s and A„. Both have the same value for the renormalized 
Coulomb repulsion /i*. For sufficiently thin films the electrons in the double 
layer travel so quickly from the superconductor to the normal metal and vice 
versa that they average over the properties of the two metals. We have essen- 
tially a new superconductor with a new averaged electron-phonon interaction A 
(in complete analogy to Cooper's arguement) 



Now we apply the strong coupling gap equation (J7J to this artificial strong 
coupling superconductor. 




- _ d s N s X s + d n N n X 
d s N s + d n N n 



'II 
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For the superconductor we chose Pb, but somewhat simplify its properties 
slightly. For the Eliashberg function we use a simple square law for Q<SId and 
express the prefactor in terms of the electron-phonon parameter A. 

We use for the Debye temperature On — 80K and for the Coulomb parameter 
ji* = 0.1. Next we calculate A (£1/) for fi; = I * 2TrksT/h = — u)j 



A(Oz) 



A 







n 2 + m 



A 1 



VL 2 



These parameters are inserted into equation JJJ) in the appendix and yield the 
dependence of T c on A. In Fig. 10 this dependence is plotted, T c as a function 
of A for (jl* = 0.1, e D = 80K and a 2 F (f2) = A (Sl/n^) 2 . The transition 
temperature of Pb, T c = 7.2K, corresponds to value of Apf, = 0.75. 



a F(Q)=x(a/n D ) 
=80K 



Project/Proximity/Supra_SC1 a 



0.4 



0.6 



0.8 



Fig. 10: The dependence of T c on the electron-phonon 
parameter A using the (linear) strong coupling gap equation. 



From Fig. 10 one obtains dT c /d\. At T s — 7.2K this slope has the value 
dT c /dX — 18.8K. Within the thin film approximation for A one finds for the 
initial slope of T c in a thin double layer: 
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dT c 




d s N s \ s + d n N n X 
d s N s + d n N n 



) 



oi- 



iis dT c 



2. 6 (X s - A„) — y 



T s dd., 



If we take A. 



0.1 and use A, 



'5 



0.75 then we find 



d s dT c 




The strong coupling treatment (in the thin film limit) reduces the prefactor 
in the intial slope to 1.7. This is a very satisfactory result. The strong coupling 
theory yields the correct intial slope. On the other hand it is also disappointing. 
It means that the superconducting proximity effect with superconductors such 
as Pb one can not be solved with the weak coupling theory. Furthermore it 
means that the simulation developed here (as all other theoretical treatments 
the author is aware of) can only be applied to weak coupling superconductors 
such as Al. But I beleave that this disappointment is more than compensated 
by the knowledge that the use of strong coupling superconductors gives an 
additional window to the study of strong coupling effects. 
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5 Conclusion 



This paper derives the transition temperature of a double or multi layer of a 
superconductor and a normal conductor numerically. The equivalence in the 
propagation of the superconducting pair amplitude and a single electron in 
Gorkov's linear gap equation is used. The single electrons act as messengers who 
carry the information about the superconducting gap (N S A (r') /tt) from one 
position-time (r', t' < 0) to another position-time (r, t = 0). This message which 
decays thermally with time as t\t (t) — J2\uj n \<n D ex P (~ ^ l*'D' ^ s integrated 
at (r,t = 0) over all start position-times (r',t') and, after multiplication with 
the BCS interaction V s , yields the new gap function A (r). At the transition 
temperature the procedure has to be self-consistent, i.e. the initial and final 
gap function have to be identical. The propagation of the single electrons is 
then quasi-classically simulated. The frame work of the calculation is the weak 
coupling theory of superconductivity. 

This numerical procedure to calculate the transition temperature of double 
or multi-layers consisting of thin films of superconductors and normal conductors 
is very flexible. The following parameters can be taken from the experiment or 
fitted during the calculation: 

• mean free path of the different metals 

• transmission through the interface 

• ratio of specular reflection to diffusive scattering at the surfaces 

• fraction of diffusive scattering at the interface. 
Furthermore it is possible 

• to vary the mean free path along the thickness of the films 

• to vary the BCS interaction NV at the interface. 

The few examples which were presented in chapter III demonstrate why the 
dirty case is so much simpler than the clean one. They also show that even for 
small thickness of the normal metal the gap parameter in the superconductor 
is not quite constant. Still the initial slope for an SN double layer follows the 
prediction of the Cooper limit. 

An important outcome of the numerical simulation is result that the nor- 
malized initial slope of an SN double layer as a function of d n at d n — does 
not depend on 

• the mean free path of the two metal 

• the thickness of the superconductor 

• a (not to large) barrier between the two metal. 
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This slope is essentially given by 



d s 



dd n 



Nn 



For the extreme weak coupling superconductor the value of T sn is 1/ (NV) . 
If one applies the numeric procedure to double layers with Pb as the supercon- 
ducting component then one obtains T sn « 4.6. This is in strong disagreement 
with the results of the few experiments which allow the evaluation of the initial 
slope. Their values for T sn lie in the range of 1.5 — 2.0. The author believes 
that the discrepancy is due to the strong coupling properties of the Pb. An 
analysis of the strong coupling gap equation in the thin-film limit confirms this 
supposition. It yields the value of about r s „ = 1.6 for a (simplified) Pb film in 
contact with a normal metal which is characterized by A = 0.1 and [i* =0.1. 
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6 Appendix 



6.1 The kernel in the clean limit 

In the clean limit the thermal Green function has the form 

G w (r, r') = - ™ exp Uk F |r - r'| ^ - M |r - r' 
27rn/ |r — r'| \ \uo\ v F 

That yields 

H u (v,v') = k B TG UJ (r,r')G: (r,r') 

or 

. . 2irk B T 1 / 2U„ 

since (r,r') depends only R = |r — r'| (using the BCS-density of state N = 
m 2 v F / (2ir 2 h 3 )). 

Without the damping the number of electrons between the radius R and 
R + dR is 

2 4^N d R= 2 -^Ndt 
nv F fr 

using dR — v F dt' . This means that (r, r') A (r') d 3 r' corresponds to an in- 
jection of 

d Z = 2 ^M1 NA (r') d 3 r'dt 

h 

electrons in the volume d 3 r' during the time dt' at the position r'. dZ is indeed a 
(dimcnsionless) number. The exponential decay exp (—2 R/v f ) corresponds 
to a decay with time since R — v F t: 

exp (-2 \uj\ R/v F ) = exp (-2 \u>\ t) 

The density of an electron at the position r and the time t — that was 
injected at (r',t' < 0) and propagates with Fermi velocity vp can described 
by the propagation density p (v F ; r,0; r', £'). Therefore (r, r') can be written 

as 

H u (r, r') = —^N (r') / dt' p (v F ;r,0; r', t') exp (-2 |«| |t'|) 



This yields the gap equation using t\t (t 1 ) from equation @ 

A (r) = V (r) / dVTV (r') / — p r,0; r', t 1 ) m (f) A (r') 

This result applies is not restricted to the clean case but applies to arbitrary 
mean free path. 
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6.2 The numerical procedure 



As shown in Fig. 2 the metal films are divided in sheets of thickness A„. Further- 
more the time developement is performed in diffusion steps of = 2X s /vp s , 
t' = mT d . Then the self-consistent gap equation takes the form 



V (z ) Td 
m=0 T ' 



(For the zero term in the time summation only half the value is taken). In the 
following we denote A (z u ) , V (z u ) , N (z v ) as A„, V V ,N V . 

6.2.1 The value of the BCS interaction {NV) S 

For the superconductor with the transition temperature T s , the density of states 
N s and the Debye temperature Q s the implicite equation 



Wrl - 1 [I/Tr] 



dt 

(NV) S ~ J ^ 
is used. 

6.2.2 Initial conditions 

At the time t = a simple gap function A„ in the superconducting film(s) is 
chosen, for example A„* = ksT s for the superconducting film(s). At the time 
t = or m = we define an occupation O v < (m = 0) of the different cells 

O v (0) = A V X V N V 

This occupation is equally divided in left and right moving electrons O v (0) 
and Ol (0) with Ol (0) = ol (0) = O v (0) /2. In the following sub-sections the 
recipe is given how to calculate from the occupation O u (m) , 0„ (m) at the 
time t = mTd the occupation O v (m + 1) , O v (m + 1). The total occupation is 
0„ (m) = O v (m) + (m). With this time devoloping occupation the new gap 
function becomes 

y oo 

A„ = — V" ?7t (mT d ) O v (m) 
K tt • £ — ' 

m=0 

This iterated gap function has two defects: (i) its shape generally does not agree 
with the original gap function A„, and (ii) the ratio of the average amplitudes 

r = ^A„^} /(Ay) will not be one. By determining numerically dr/dT from 

two iterations with the same initial gap function and two temperatures T and 
T + Ta the temperature is adjusted, using Newton's extrapolation method. 
After a few iterations ( A„) becomes sufficiently close to (A„) and the adjusted 
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temperature is the transition temperature of the multi-layer. The iteration is 
completed when 

v Zs v < 10- 5 
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6.3 Diffusive and ballistic propagation 

The important task is to devise a simple fast procedure that describes the bal- 
listic propagation of the electrons for distances shorter than the mean free path 
I and the diffusive propagation for distances larger than /. It helps consider- 
ably that only the propagation in z direction has to be modeled properly (as 
long as no magnetic field perpendicular to the film is applied) . We consider the 
electrons in a thin layer of thickness dz in the interval (z, z + dz). Half of the 
electrons have a positive z component v z — vf cos 8 of the velocity. As long as 
they are not scattered their average velocity in z direction is 

C^ 2 2it sin Ov f cos Odd l 
J* /2 2irsm0d0 2 

We take this as the minimum requirement for the ballistic simulation. 

The simulation of the diffusion in z direction is rather straight forward. At 
the time t = 0we have the initial occupation 0„ (0). 

Let us first consider the diffusion in one dimension. Here the electrons have 
cither the velocity +vp or —vp. The size of the cells is A and an electron needs 
the time eo—X/vp to cross a cell. We divide the initial occupation O v (0) into 
O v (0) = O v (0) = O v (0) /2 for the left and right moving electrons. When 
the electrons reach the boundary of the cell they will be partially transmitted 
through the boundary with the probability p and partially reflected with the 
probability (I — p). This yields the rule how of one obtains from the occupations 
at the time t — mso the occupation at the next time step t = (to + 1) £q 



O v (to + 1) = pOv-x (to) + (1 - p) O v (to) 



O v (m + 1) = pO v+1 (to) + (1 - p) O v (to) 
This yields a one-dimensional diffusion with the diffusion constant D = 

1 p A 2 

2 1-p e ' 

Ballistic propagation requires setting p almost equal to 1 . In this case almost 
all the O v (to) electrons move from cell v to cell (v + 1) during the time £o- 
This means that they propagate the average distance A = vfEo during the 
time £o- Therefore this model does not fulfill the basic requirement for ballistic 
propagation in three dimension that (v z ) = \vp. 

A three-dimensional diffusion can be obtained by a sequential propagation 
in x, y and z direction, each for a time of £o with the velocity vp. This yields a 
diffusion constant D = 5 ^ and triples the average time for the diffusion in 
z direction. Since the electrons propagate only during every third of the interval 
3e in z direction they propagate the distance A during the time 3e , i- e - their 
average velocity in z direction is only (v z ) = v F /3. 

We can simulate the average diffusive and ballistic propagation of the elec- 
trons in z direction by propagating every other time interval £o in z direction. 
Then the time step is Td — 2eo- In this case the diffusion constant is D = 
and the ballistic propagation yields (v z ) = vp/2 as required. 
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It should be mentioned that it is essential that the electron density is divided 
into (at least) two components, one for motion in the +z and the other for 
-z direction. A single density component with hopping to neighbor places 
yields only small diffusion constants of D = and can't describe the ballistic 
propagation at all. 

For the normal conductor the same time element Td is used to simulate 
the propagation. The thickness d n is divided in cells (or layers) of thickness 
An = VF, n £o = VF, n Td/2. This synchornizes the diffusion in the whole double 
layer. 

The transparency p of the cell walls is obtained form the experimental con- 
ductivity a of the films, where a m = 2e 2 N m D m or 

p. &m Dm 

where m stands for s or n. 
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6.4 Interface between two films 



The transmission of electrons through an interface between two metals (which 
we denote with S and N) is only in exceptional cases equal to 1 . If for example 
the Fermi wave number kp iS is larger then kp^ n then any electron in S whose 
component k p parallel to the surface is larger than kp^ n cannot cross the interface 
because afterwards it would have an energy of a least (hk p ) 2 /2m which is larger 
than the Fermi energy Ep n — {hkp ^ j Ira in the normal conductor. An 
electron in N with Fermi energy would not violate the conservation of energy 
when crossing the interface. However, a plane wave which crosses a step in the 
potential energy is partially reflected. Therefore the transition probability is 
less than 1 for any electron. If one averages the transition probability of all 
these electrons (to cross the interface from N to S) one finds 

Tn->s = f i^r^- - > where 

4 (V( X + 1 )Y (x + 6) - (V^) 5 - 10* - 6 
I W = 15 " ? 

For small x the asymptotic expansion is / (x) ^ (l — *>/%)■ 

The detailed balance requires that in equilibrium the number of electrons 
which cross from S to N is equal to the number of electron which cross from 
N to S. Let us assume that the electron distribution is in equilibrium and we 
consider an interface S/N. O s (m) and O n (m) are the occupations in the cells 
on the left and right side of the interface. The transmission coefficients are by 
T sn and T ns . Then the occupation at the time (m + 1) is 

0~ s (m + 1) = T ns <T n (m) + (1 - T sn ) ol (to) 
(m + 1) = T sn (f s (to) + (1 - T ns ) On (to) 

In equilibrium one has O s , n = O s , n = -O s n and O s „ (m + 1) = O s n (m). This 
yields 

T ns On (TO) = T sn O s (to) 

Since O s _ n — X s ^ n N s _ n one obtains finally 

T ns X S N S 

If one considers real metals a considerably more difficult situation arises 
when the superconductor has a mass enhancement of the density of states (as 
most superconductor have, in particular the strong coupling ones). We return 
to the mass enhancement below. However, independent how complicated the 
individual transmission probabilities are, the detailed balance will always apply. 
In our simulation we use T ns < 1 as a fit parameter and calculate T sn using the 
detailed balance. 
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6.5 Strong coupling effects 

In appendix (6.1) the kernel H u (r, r')was derived for free electrons 
H„(r,r') = k B TG u (r,r')G* w (r,r') 

where 

Gui (r, r') - - 7rl ^ jr exp (ik F |r - r'| ^ - M |r - r>\) 

2Tth z \Y — r'| \ |w| / 

G w (r, r') is a Fourier transform of G w (k) 

G w (k) = — ^~ 

- £fe 

In the renormalized strong coupling case one has 

1 1 



G w (k) = 



1 + A iuj — £fc 



where (1 + A) is the "mass enhancement" of the electrons at the Fermi surface 
due to the electron-phonon interaction. Here A is defined as 



A = 2 f D a 2 FH — 
Jo u 



and a 2 F (ui) is the Eliashberg function of the electron-phonon interaction. In 
performing the Fourier transform one obtains 



G " (i?) =2^ 6XP 



UJ \(jj\ 

ik F R— - -jr-i-R 

M nv* F 



(using |r — r'| = R). Compared with the free electron Green function the Fermi 
velocity is now renormalized. 

This yields for the function (R) 

27rk B T N 1 / 2H 

using the bare BCS-density of state N = m 2 VF /2ir 2 h 3 and the bare Fermi 
velocity Vf ■ 

Now we use the same argument as before: The term (r, r') A (r') d 3 r' 
corresponds to an injection of 

dZ = -±—^L N A (r') d 3 r'dt 

(l + A) h y J 

electrons in the volume d 3 r' during the time dt at the position r' which propagate 
with v F . The factor 1/ (l + A) arrives from the ratio of v f /vf- The exponential 
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decay exp (—2 \ui\ R/v* F ) corresponds to a decay with time as exp (—2 |w| R/v* F ) = 
cxp (—2 \oj\ t 1 ) since R — v* F t' . Now we replace the BCS interaction V by an 
effective interaction strength (A — fj,*) /N. (a 2 F (oj) and therefore A contains 
the bare density of states iV as a factor) 
The gap equation is then 

a(r) _ iM^fi^ f\, lMlTm) 

A test for constant A using J d 3 r' p {v* F \ r,0; r', t') — 1 yields 

A = ~ h TTa J_ 00 VTc ( t ^ dtA 

This yields McMillan's first approximation for the T c of a strong coupling su- 
perconductor 

(1 + A) = "f?" _J_ 
6.5.1 Strong coupling gap equations 

For the strong coupling superconductor one has to use the Eliashberg theory. It 
replaces the "one" gap equation by series of gap equations at different Matsubara 
frequencies. Equation Q 42 defines the gap A (u>i) at uji as a function of the 
gap at all other Matsubara frequencies u)j. 

A (a*) = 2 -^f- J2 {A [uh - u>i) - A 4 *} ^A (7) 
7rfc B T ^ ui t . 

^ = + —f— m A - w o 

no/" 00 , wa 2 f(w) 
A (Wj - WjJ = 2 / dw— 

JO OJ + [U>i — OJj) 

The effective BCS interaction NV is replaced by {A (a;, — a/j) — //*}. The 
electron-phonon parameter A (wj — ojj) is determined by the Eliashberg function 
a 2 F(w). The denominator 1/2 \u n \ is replaced by the dressed Matsubara fre- 
quencies 1/2 \cJj\. Here is is ulj — ujj (l + A] where A is an average over (2j + 1) 
different values of A (0^), (Oy = v2nkBT/H,0 < ^ < j). 

One can derive from the energy dependent gap equation a renormalized gap 
equation by the following simplifications 

• A (uii) => A = const for \uji\ < £Id 
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• A (wj) => for \u>i\ > Op 

• one replaces A — u)j) => A (0) = A for , < 12 £> 

Then one finds that uJj is renormalized as = Uj (1 + A) where (i + A) 
is the electron-phonon renormalization factor. Using these simplifications the 
renormalized gap equation takes the form 




as before. 
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